Fortran 浮点数的逻辑运算

! 注意浮点数的运算误差
a = sqrt(b)**2 - b
if (a == 0.0) ...
if (abs(a) <= e) ...